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ABSTRACT 



A toy model of a disk undergoing steady state accretion onto a black hole is presented. 
The disk is in a hydrostatic equilibrium for all radii r > rin, with the inner disk radius 
located between the marginally stable and marginally bound orbits: r„is > rin > r,nb- 
Matter flows from the disk through a narrow cusp at rma and falls freely into the black 
hole, carrying with it no thermal energy. At radii larger than rout the disk is assumed 
to radiate away all locally generated heat, and therefore the disk is geometrically thin for 
r > rout- We assume that no heat generated in the inner disk, with rout > r > ri„ can 
be radiated away, i.e. the disk is 100% advective, and it becomes geometrically thick in 
this range of radii. All enthalpy of the thick disk is used up to press the inner disk radius 
towards the marginally bound orbit, and to lower the efficiency of conversion of accreted 
mass into radiation generated only for r > rout, by assumption. 

Conservation laws of mass, angular momentum and energy make it possible to calcu- 
late the inner thick disk radius ri„ for any specified value of its outer radius rout- As the 
nature of disk viscosity is not known there is some freedom in choosing the shape of the 
thick disk, subject to several general conditions, which include the hydrostatic equilibrium 
everywhere for r > rin - The main purpose of this toy model is to emphasize the effect the 
disk thickness has on lowering the energetic efficiency of a black hole accretion. 
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1 Introduction 

Following the publication of the first paper on the advection dominated ac- 
cretion flows by Narayan and Yi (1994) over one hundred papers appeared 
on this subject, with ever more sophisticated physics and ever more sophisti- 
cated treatment of geometry (cf. Gammie and Popham 1998, and references 
therein). However, there is no fully two dimensional treatment |^ published 
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so far, and there is no clear link between the recent papers and those writ- 
ten about two decades ago on the theory of thick disks (e.g. Jaroszynski, 
Abramowicz and Paczynski 1980, Paczyhski and Wiita 1980, and references 
therein). The purpose of this paper is to present a toy model of a disk 
accreting onto a black hole but not radiating, i.e. advection dominated, 
presented in the spirit of the early 1980s, which seems to be simpler and 
more transparent than the spirit dominating the late 1990s. 



2 Thin Pseudo-Newtonian Disk 

I shall use pseudo-Newtonian gravitational potential (Paczynski and Wiita 
1980): 

GM , _d-^ _ GM _ 2GM 

E^g' dR' {R-Rgf ^5-^2-' (1) 

which has the property that a test particle has a marginally stable orbit 

at Rms = 3Rg, and marginally bound orbit at Rmb = ^i?^, just as it is 

the case with particles orbiting Schwarzschild black hole. R is the spherical 

radius. Throughout this paper cylindrical coordinates {r,z) will be used, 
with/2= (r2 + z2)i/2_ 

In the following we consider a thin accretion disk fully supported against 
gravity by the centrifugal acceleration. The disk is in the equatorial plane 
of the coordinate system. It follows the motion of test particles on circu- 
lar orbits, with the rotational velocity v{r), angular velocity 0(r), specific 
angular momentum j(r), and total specific energy e(r) given as 
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where we adopted rg = Rg. It follows that 

dr dr ^ ^ 

The inner edge of a thin accretion disk is at 

fin ~ fms ~ ^fgi (4) 

where the binding energy and the specific angular momentum reach their 
minima: 

Cms — ~Y6' 

1/n GM GM 
jms = 1.5 X 6^/2 X ^ 3.674 . 6 

c c 

The matter falls freely into the black hole once it crossed the r^s , conserving 
its angular momentum and total energy. 

The total thin disk luminosity is given with the formula: 

2 

Ld = Mcin = Mems = (-m) ^, (thin disk), (7) 

where I adopt a convention that M < for the accretion flow. 

In a steady state accretion of a thin disk the equations of mass and 
angular momentum conservation may be written as 

M = const. (8a) 

j = Mj + g = const. (86) 

where g is the torque acting between two adjacent rings in the disk. With 
no torque at the inner edge of the disk at Vin = Vms we have 

g = (-M) (j - Mn) = (-M) (j - jms) . (9) 

The rate of energy flow in the disk is given as 

E = Me + gn. (10) 

This is not constant, as the disk must radiate energy dissipated by viscous 
stresses in order to remain thin. The energy balance equation may be written 
as 

dE / dn\ / ■ \ , , / dn\ , , 

. 2„ = __ = , = (_M) - ,-.„) (--) . (11) 
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where F is the energy radiated per unit disk area from each of its two 
surfaces. The eq. (11) may be integrated from rms to infinity to obtain the 
same result as that given with the eq. (7). 

Note, that disk viscosity was not specified. All we needed were the 
conservation laws and the assumption of a steady state accretion of a geo- 
metrically thin disk. Small geometrical thickness implied that no significant 
internal energy could be stored within the disk. 



3 Thick Pseudo-Newtonian Disk 

1 shall consider now geometrically thick, axisymmetric disk, with the surface 
given with the relation Zs{r), where (r, z) are the two cylindrical coordinates, 
and we have 

R=(r^ + z^y^\ (12) 

The disk is assumed to be in a hydrostatic equilibrium. This implies that 

the three vectors have to balance at its surface: gravitational acceleration, 
centrifugal acceleration, and the pressure gradient divided by gas density. 
This implies that angular velocity fig at the disk surface is given as (cf. 
Jaroszyriski et al. 1980, Paczynski and Wiita 1980): 

which is equivalent to 



R<,J \ r dr 



where subscript 's' indicates that the corresponding quantities are defined 
at the disk surface. 

Thick disk has a large radial pressure gradient, and it remains in a hy- 
drostatic equilibrium for r > r^n, and it has a cusp at rin- The matter 
falls freely into the black hole inwards of rin, which is located between the 
marginally stable and marginally bound orbits, i.e. 

'^mb ^ fin "^ms (-^'^) 

where a no torque inner boundary condition is applied. I shall not analyze 
the transition from a sub-sonic accretion flow at r > Tin to a supersonic flow 
at r < Tin (cf. Loska 1982). 
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The total disk luminosity is given with the formula similar to eq. (7): 
La = Me,n = (-M) GM{nn-1r ) ^ ^^^.^j^ ^^^^ 

^\Jin fg) 

The specific binding energy at ri„ has the range 

0> ein> -— for rmb < nn < rms- (16) 

Note, that while the total luminosity of a thick disk may be larger than that 
of a thin disk, the energy radiated per each gram of matter accreted into 
the black hole is smaller in the thick disk case. 

4 Thick Advection Dominated Disk 

A thin accretion disk has to radiate in order to remain thin. If the disk is 
not radiating away the energy which is generated by viscous stresses then 
it must be thick in order to accommodate this energy in its interior. Let us 
consider a disk which is thin for r > Tout-, but is thick for Tout > f > Tin- 
Note, that the disk is also thin at its cusp, i.e. at r = rj„. Our new disk 
radiates energy only where it is thin, i.e. for r > rout, where the disk surface 
brightness F is given with the eq. (11), and the total disk luminosity is 
given by the formula similar to the eq. (7) and (15): 

La = Mcin. (17) 

All this luminosity is radiated by the thin disk at r > routt and none is 
radiated by the thick disk at ro^t > r > rin, by assumption. 

The structure of the thick disk, extending from rin to rout is of interest 
for us, as by assumption the disk is 100% advection dominated in this range 
of radii, i.e. the mass, angular momentum, and energy entering it at rout 



must come out of it at rin- 

Min = Mont- (18a) 
Jin = (Mj +g), = (Mj +g) = Jout, (186) 

V / in V / out 

Bin = (Me + gn) = (Ue + gO) = Eout (18c) 

\ / m \ / out 

with j and O having their 'Keplerian' values (cf. eqs. 2) at the two ends, as 
the disk is thin at both ends, and there is no torque at rj„, i.e. gin = 0. 
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Figure 1: The variation of the inner thick disk radius rin and the corre- 
sponding binding energy is shown as a function of the outer thick disk 
radius rout- fg = 2GM/c^ is the black hole radius, and e^s = — c^/16 is 
binding energy of a test particle with a unit mass orbiting the black hole at 
the marginally stable orbit with radius Vms = ^fg- 

The torque at rout follows from the conservation of angular momentum 
(18b): 

gout = {-^) {jout - Jin) , (19) 

which, together with the energy conservation law give us: 

(^out (^in ~ {jout jin) ^out- {^^} 

Note, that all quantities in the eq. (20) arc unique functions of cither rj„ or 
'''out (cf- eqs. 2), and therefore the eq. (20) gives the relation between r^ 
and rout- If the outer radius of the advection dominated disk is specified, 
then the inner radius of the disk can be calculated with the eq. (20), and 
the eq. (17) gives the total luminosity of the thin disk which is assumed to 
extend from rout to infinity. The relation between rout and rj„ is shown in 
Figure 1. The larger is Vout, the smaller is Vin, and the less energy is radiated 
per one gram of accreted matter. 
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The variation of the inner thick disk radius rin, and the binding energy 
at the inner thick disk radius e^n with the outer thick disk radius rout is 
shown in Fig. 1. Notice, that for very large rout we have asymptotically 

!i!i~2 + 3^, for ^ » 1. (21) 

fg fout f^ms fout fg 

We also have 

eout^~ — , 3eo„t ~ — , for — » 1. (22) 

^ font ^ fout fg 

So far we used only conservation laws to constrain our thick advective 
disk. If we want to find the disk shape we must make some additional 
assumptions. There are two general inequalities which must be satisfied by 
the matter at the disk surface (cf. Jaroszyhski et al. 1980, Paczyhski and 
Wiita 1980): 

supplemented with the condition of hydrostatic equilibrium at the thick disk 
surface, as expressed with the eq. (13b), and the conditions that the disk 
must be geometrically thin at rj„ and Tout- 

There is a lot of freedom in choosing thick disk structure that satisfies 
all the conditions listed in the previous paragraph. For our toy model we 
adopted 

1.5/a 

(24) 



a / \ 3a 

l + 6( — -1) +6^-1 



Js = 3 
and 

djs f^" 1 dj^ 
es = ein+ ^s-j-dr = ein + '^'f'dr- (25) 

Thin disk conditions at rj„ are satisfied automatically with the eqs. (24) 
and (25), and we have 

{GMrfS^' _ GM{r,^-2r,) 

~ r- -r ' ~ ')(r- -r V ' ^ ' 

I m I g m I g ) 

The parameters a and h have to be adjusted so that thin disk conditions are 
satisfied at Tout-, where the eqs. (24) and (25) must give 

_ {GMrlJ^'^ GMjrout -2fg) 

Jout- ' ^out- ^, ,_r)2 ■ ^"^'1 

I out I g V ' out I g ) 
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Figure 2: A cross section of the toy disk model accreting into a black hole. 
All energy generated inwards of Vguf — 100 rg is assumed to be advected with 
the accretion flow, and none is radiated away. At large distance, r > rout: 
all energy generated by viscosity is assumed to be radiated away, hence the 
disk is geometrically thin for r > rout- The inability to radiate away energy 
for r < rout forces the disk to become thick, and pushes its inner disk radius, 
rin = 2.026 rg, close to the marginally bound orbit, r^b = 2 rg, where the 
binding energy is only |ei„| !^ 0.006 (P. 

As an example a toy disk model with rout = 100 rg is shown in Fig. 2 
and Fig. 3. The inner disk radius r^n = 2.026031... rg was obtained from 
the condition given with the eq. (20). The eqs. (27) were satisfied for 
a = 1.32048586889... and b = 0.000000442466655.... 

5 Discussion 

The shape of our toy disk as shown in Fig. 2 is an artifact of the strong 

assumptions made in this paper. In particular, a rapid transition from a thin 
disk to a thick disk at rout is a direct consequence of the assumption that 
100% of all energy dissipated locally is radiated away locally for r > rout, 
while none is radiated away for r < rout- In any realistic disk there will be 
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Figure 3: The inner parts of the toy disk model presented in Fig. 2. The 
dots on the equatorial axis indicate the location of the marginally bound 
orbit at r^f, = 2 r^, and the marginally stable orbit at r^s = 3 r^. The disk 
becomes geometrically thin at rj„ = 2.026 r^, and cold matter falls freely 
into the black hole at r < r^n- No thermal energy is advected into the black 
hole in this model. 

partial radiation at all radii, and the variation of the efficiency is likely to 
change gradually, with no abrupt changes in the disk thickness. 

While the detailed shape of the thick disk must be uncertain as long as 
we have no quantitative understanding of disk viscosity, the formation of a 
thick disk and pushing its inner radius towards the marginally bound orbit 
is a very general property. It was noticed two decades ago with supercritical 
accretion disks of Jaroszyhski ct al. (1980) and Paczyhski and Wiita (1980). 
In our toy model the inability of the disk to radiate energy dissipated within 
font > r > Tin forces the disk to become thick, and pushes its inner cusp 
towards Vmb^ lowering the efficiency with which rest mass is converted into 
radiation. The farther out the Tout is assumed to be the closer ri„ is pushed 
towards the r„ib hi order to lower the efficiency. 

While we retain the term 'advcctive' to describe our thick disk, it should 
be stressed that as we require the disk to be thin at its inner cusp at rj„, there 
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Figure 4: The variation of angular momentum with radius for the surface of 
the thick disk presented in Figures 2 and 3 is shown with the sohd hne (js). 
The variation of 'Keplerian' angular momentum with radius is shown with 
the dashed line {Jk)- Note that thick disk angular momentum is almost 
constant for r < 10 r^. 

is no advection of any heat into the black hole, as all enthalpy of the thick 
disk has been used up to press its inner radius towards the rmb- However, the 
kinetic energy of thick disk matter at rin is much larger than is the kinetic 
energy of a thin disk, which must have its inner radius located at rms- In 
the pseudo-Newtonian potential the ratio of the two is v^^/v^g = 8/3 (cf. 
eq. 2a). 

There is a somewhat paradoxical aspect of our toy model. Some kind 

of viscosity has to reduce angular momentum from jout at the outer disk 
boundary to jin at its inner boundary, and this must generate either heat 
or some other form of internal energy which puffs up the disk, and it would 
seem that this energy has to be advected through rj„ into the black hole. 
Figure 4 shows the distribution of angular momentum for our thick disk 
model (solid line), and for the 'Keplerian' angular momentum given with 
the eq. (2c) (dashed line). It is clear that thick disk angular momentum is 
almost constant for r < 10 r^, and that all dissipation takes place only in 
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the outer parts of the thick disk. All the same the entropy must be higher 
at Tin than it is at rout- However, high entropy does not have to imply high 
internal energy if the gas density is low at Tin- 

Obviously, our requirement for the disk to be thin at rj„ is ad hoc, mo- 
tivated by our desire to have as simple structure as possible. Our toy model 
may require uninterestingly low accretion rate to satisfy all the conditions 
that are imposed at its inner boundary: hydrostatic equilibrium, low geo- 
metrical thickness, high entropy and low internal energy. It is very likely 
that at the accretion rate of any interest some of these conditions are broken. 
It is likely that the cusp at rin opens up considerably, the speed of sound is 
not negligible compared to rotational velocity, and the transonic flow car- 
ries a non trivial amount of internal energy into the black hole (Loska 1982). 
However, the conservation laws do not require the advected thermal energy 
to be large, as demonstrated by our toy model. 

A generic feature of any thick disk is the formation of a narrow funnel 
along the rotation axis. It is far from clear how realistic is the presence 
of the funnel, as some instability might fill it in, and make the accretion 
nearly spherical near the black hole. On the other hand, if the funnel forms 
and lasts, it may help collimate a jet-like outflow. Unfortunately, our lack 
of quantitative understanding of viscous processes in accretion flows makes 
it impossible to prove what topology of solutions is realistic under which 
conditions. The main virtue of the toy model presented in this paper is 
the set of assumptions that was made; this set is very different from the 
assumptions which are used in the recently booming industry of advection 
dominated accretion flows. While the assumptions adopted in this paper 
are ad hoc, so are the assumptions adopted in any thick disk models. In 
particular, there is nothing 'natural' about the popular assumption that 
the accretion flow is self-similar. Another popular assumption: constant 
'alpha' parameter is ad hoc as well. It is useful to read the old paper by 
Galeev, Rossncr and Vaiana (1979) to realize that the very concept of the 
'alpha' parameter is ad hoc. Disk properties which depend on any ad hoc 
assumptions should not be taken seriously. 
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